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Abstract 

In this Letter we propose a scheme to generate and detect entanglement 
between charge states in superconducting nanocircuits. We discuss differ- 
ent procedures to discriminate such entanglement from classical correlations. 
The case of maximally entangled states of two and three coupled Josephson 
junctions is discussed as example. 
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The phenomenon of entanglement, probably one of the most striking feature of quantum 
mechanics Q, appearing as a consequence of the superposition principle in the presence of 
composite systems, is the main ingredient in all known examples of quantum speed-up in 
quantum computation and communication ||2|||. This has prompted intense experimental 
efforts towards its generation and detection, most notably with photons ||, cavity QED 
systems ||, ion traps @, and coupled quantum dots @. Although in condensed matter it 
is common to encounter correlated many-body states, it is difficult to isolate the different 
subsystems while maintaining their entanglement. This problem is recently attracting a lot 
of attention and several solid state devices have been suggested. They are based on the 
phenomenon of Andreev reflection in hybrid normal-superconducting systems |||9j or on the 
coupling of mesoscopic Josephson junctions with superconducting resonators [|K],|ll[ . 



In this Letter we propose an explicit experimental scheme to generate and detect entan- 
glement in superconducting nanocircuits. Astonishing progresses have recently been made 
in the control of the coherent evolution of such systems ||T2"| , |TB|| , which have been proposed as 



promising candidate to realize a quantum computer [p~4| — pLTfl . In our proposal we specifically 
address two important issues, i) We describe how to measure the entangled states once the 
subsystems have been separated and what is the effect of a possible residual coupling on the 
outcome of the measurements, ii) We show how to distinguish between entangled states and 
statistical mixtures or product states. Our setup is based on a modification of the device 
used by Nakamura et al. |T3|JT8| and, we believe, it is amenable of experimental verification 
with present days technology. 

As sketched in Fig.|I], we consider two superconducting Single Electron Tunneling (SET) 
transistors (labelled by a, b) coupled by a small, Josephson junction. By choosing appro- 



priately the working point of the device [13| (see below), coherent Cooper-pair tunneling 



takes place only across the left and the coupling junctions, while quasi-particle tunneling 
is important across the right junctions. With these observations the Hamiltonian can be 
written as a sum of three contributions H = Hq + H qp + Ht- Its coherent part, Hq, is given 
by 



H = [ E ch(ni - n X)i ) 2 - Ej cos 2ip { 



i=a,b 



+ Ei ab \n a - n x , a )(n b - n x , h ) + cos2(<p a - <p b ) (1) 

Here, E c h is the charging energy, n Xti are offset charges induced by external voltages and Ej 
is associated to the Josephson tunneling. The phases ipi and the number of charges on the 
islands rii are conjugate variables [ipi,rij] = i 5ij. The state of each individual SET can be 
manipulated by varying n x with a suitable choice of voltage pulses, for example by putting 
the left junctions at resonance for Cooper-pair tunneling. However, in order to produce 
entanglement, we must be able to perform transformations other than local ones. This 
requires a controllable coupling between the SETs. In our proposal such coupling, Ej (<&), 
is provided by a SQUID, pierced by an external flux $. The SQUID capacitance, which 
we assume to be much smaller than all the other capacitances in the setup, gives rise to an 
electrostatic coupling, E^. The term H qp = EaEfccr £ kalkaalkaa ( a = R a ,Rb,h,h labels 
the right electrodes and both islands respectively) describes the quasi-particles. Here, T k(Ta 



ilkaa) creates (destroyes) a quasiparticle with momentum k and energy e ka = y£f Q + A 2 (£ k 
is the single-particle dispersion, A is the superconducting gap and o the spin label). Finally 
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where T kq is 



the tunneling Hamiltonian H T is H T = Ei=a,& [e l<ft Y,k q aTk q ll (TRi 7 q ai l + h.c 
the amplitude for quasi-particle tunneling. By fixing the transport voltage \eV\ ~ 2 A + E ch 
and due to the Coulomb blockade E ch ^> Ej, it suffices to consider only three charge states, 
n = 0, 1,2 and two quasi-particle tunneling rates, T^ b and T^ b . The latter are related to 
transitions across the right junctions with the charge on the island changing as n = 2 — > 
n — 1 and n = 1 — > n = respectively (in the regime we consider, ~ r^ a//b ^ ~ T ) [I9|. 

The other tunneling rates are exponentially small. We furthermore assume T -C Ej. 

The dynamics of the system can be described by a master equation for the density matrix 
representing the charge state of the two islands |19| 



dtP = i[p,H Q ]+± £ £ r«{2L«p(4 l) ) t - {Lf)'Lfp-p{Lf)'Lf) (2) 

Z i=a,bn=l,2 1 ' J 



where L$ — \n — l)(n|$ is the Lindblad operator corresponding to the quantum jump 
n — > n — 1 for the i-th island. The states |0) and |2) of the two SETs are involved in 
generation of entanglement, while quasi-particle transitions through the states |1) allow to 
perform the quantum measurement. When such devices are employed as qubits states |0) 
and 1 2) are the computational states. 

Ideally the detection of entanglement goes through the following steps, i) Prepare the 
entangled state by means of manipulation of the gate voltages; ii) Switch off the coupling 
between the two qubits; iii) Perform the measurement. 

Preparation - We illustrate one of the possible procedures to prepare the singlet state, 
!*(-)} = (|02) - |20))/>/2. The two junctions are initially kept off de generacy with the 
initial state given by 1 02) . Then, by switching on the Josephson coupling Ej 1 ^, and slightly 
shifting the working point of the two qubits such that n x>a — n Xj b — Ej^ /E C h, for a time 
St = ir/(V2E { j ab) ), the desired singlet state is obtained. In a similar fashion it is possible to 
generate other maximally entangled states. After the preparation, the coupling between the 
two SETs should be switched off. 

Equally important to the ability to generate entanglement is the possibility to detect it. 
Below we discuss two possible detection schemes. 

Method 1 - The success in the preparation of against decoherence can be tested, 

as proposed in Ref. ]2U|, by measuring the correlation between the integrated quasiparticle 



current signals of the two SETs, I a /b = Jo™ I a /b(t) dt, with the integration time, t m , chosen 
to be much longer than The correlator is given by (e is the electron charge) 

C = e- 2 {IJ b ) = e~ 2 / dti / dt 2 {Ia{ti)h{t 2 )) (3) 
Jo Jo 

where the two-time average {I a {ti)h{t2)) can be obtained with the help of the quantum 
regression theorem applied to eq. (0). Following the procedure of Nakamura et al. [13| , 
the experiment can be performed by means of repeated preparations of the initial state and 
measurements of the correlator. A value of C different from zero indicates that current is 
flowing in both channels. This cannot happen if the system is prepared in the state: 
if island a relaxes from |2) to |0), then b does not and viceversa. 

Complete anticorrelation between the two currents is not enough to come to the con- 
clusion that the charge state of the two islands is entangled. A statistical mixture of 
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the form p mix = |(|02)(02| + 1 20) (20 1 ) or any convex combination of the kind p p = 
PPmix + (1 ~ ]?)|\I/ < '~-')(\I / (~- ) | would also lead to the same result for C. In order to discrimi- 
nate between quantum vs. classical correlations we can proceed in the spirit of experiments 
testing Bell's inequality [l|]. The basic idea is that if the system is in state the result 

C = is obtained also if, before the detection stage, one carries out a (further) bi-local 
unitary operation. For example, one can rotate the state of the two qubits by the same 
angle 9 bringing both islands to degeneracy for a time interval At = 29/ Ej. This leaves 
|\l/( _ )) unaffected. On the other hand when the system is in a state like p m i X bi-local unitary 
transformations will in general produce a non-zero population in 1 22) therefore giving rise 
to a non-zero C. This is the key ingredient of our scheme to discriminate entangled states 
from classical mixture. 

In the ideal case we can think to switch off the Josephson interaction energy during the 
measurement, thus "separating" the two subsystems. If this can be done and if, in addition, 
also the Josephson energies of the two transistors can be set to zero, then the master equation 
can be solved exactly. Starting with the state p p , the integrated current correlation is given 
by 

Tt 



C p = psm 2 (29) 



I _ e -rt m 1 b rn ^-rt r . 



(4) 
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The two cases of the singlet state and of the complete mixture are obtained for p = and 
1, respectively. Note that for p = 1 an oscillatory behaviour is obtained, whose visibility is 
reduced for p < 1. Furthermore only when p = the correlator C vanishes. 

In order to ascertain the validity of this proposal it is important to check how the residual 
couplings affect the measurement. Experimentally a non-zero electrostatic coupling will be 
always present; however, as long as it is small compared to E c h, it has no effect on the 
measurement. More important is the role of Ej and of the coupling Ej 1 ^; nevertheless, 
even if residual Josephson coupling energies are present during the measurement, the singlet 
and the ] statistical mixture give rise to very different correlation signals. In Figs.|2|,0 the 
correlation C is shown for t m = 3r _1 , for both the p = and p = 1 cases, and for different 
values of the charging energy of the SETs and E^j b \ respectively. During the measurement, 
the Cooper pair states n = 0, 2 should be kept off degeneracy. The amplitude of the 
oscillations in the mixed case depends on the charging energy E c h (for finite Ej) because 
the coherent oscillations in the two SETs may significantly change the state during the 
measurement. Differences between the pure and mixed state in Fig.^| can be detected down 
to values of E c h/Ej > 1 (see the dashed line in Fig.|]). For lower values of E c h/Ej, it is 
impossible to discriminate between the p = and p = 1 cases because of the extra Cooper 
pair tunneling. It is theoretically desiderable and possible experimentally to fix the ratio 
E c h/Ej 3> 1, therefore this is not a limitation for the feasibility of our proposal. The signal 
does not show substantial variations with increasing the Josephson coupling E^ either (see 
Fig.|3]). The singlet state always gives rise to a correlation function which does not change 
with 9. Additional errors, deriving from not perfect gating, are possible; but they are not 
specific for the measurement of entanglement and it was aready verified that they do 
not affect the preparation of the state. One can conclude that our proposal is robust against 
various non-ideal situations that can be encountered in an experiment. 

In order to perform current correlation measurements one needs detectors sensitive to 
single electrons, which is technologically very demanding. One possibility could be to amplify 
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the signals coming from the two qubits before correlating them. A similar procedure has been 
recently employed to measure the correlation between the two output ports of a fermionic 
beam splitter |2T] . 



Method 2 - The protocol described so far is particularly suitable to verify that a given 
maximally entangled state has been created and to detect if decoherence mechanisms have 
spoiled its coherence. A different approach, which does not require the technologically 
demanding current correlation measurement (nor an a priori knowledge of which state has 
to be detected) can be designed as well. The basic requirement of this approach is the ability 
to implement a two-qubit gate, E2J . 

To this end we consider the following sequence of voltage pulses, i) Bring the 6-island to 
degeneracy with a voltage pulse of duration ti = n/2Ej; ii) Set n~\ = n~\ = 1 + E^ /2E c h 
and let the system evolve for a time t\ = tt/AE^; iii) rotate again the 6-island of an angle 
37r/4 keeping it at degeneracy for a time £3 = 3n/2Ej. In the computational basis this 
corresponds to the operation 



(5) 



which is a C-Not gate apart from a phase factor. It is easy to verify that the deviations 
from ideal gate operation due to the presence of the electrostatic coupling during single 
qubit operations are 0{E < f b ' > / Ej) 2 . The gate Q can be used to transform entanglement into 
coherence of the control qubit, thus providing a mechanism to discriminate entanglement 
from classical correlations. Indeed we have: 

g\&-))(&-)\g — > I(|0> - i|2»«0| + i(2\) <g> |2)(2| (6) 

GpmixG — - -I® |2)(2| (7) 

Note that the target qubit (b) always factors, while a is left in either a coherent or a 
completely incoherent superposition of the basis states, depending on the initial presence of 
entanglement. At this point, one can reveal the coherence of the state of the control qubit 
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by performing a current measurement as in [13|. If the excess current shows oscillations with 



respect to a varying rotation angle, then the state is coherent and therefore the initial state 
was a singlet. On the other hand, if the excess current coming from island a is constant one 
can infer that the initial state was not entangled. 

We mention that a Bell state analyzer can be easily implemented with this setup. Indeed 
by applying Q followed by a further 7r/4 rotation, the four maximally entangled Bell states, 
|$(±)) = ^-(|00) ± |22)) and |^ (±) ) = ^(|02) ± |20)), are transformed into the basis state 
1 20), 1 00), 1 02) and |22), respectively, so that they can be distinguished by two local single- 
qubit current (or charge) measurements. 

All the preceding discussion could have been equally well phrased in terms of charge mea- 
surements instead of current ones. In this case, the right junctions are no longer necessary, 
while the remaining single Cooper pair boxes should be coupled to two electrometers |23| 



which, in turn, induce a dephasing and a mixing among the qubit states EM. The procedure 
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can be easily generalized to include more qubits. With a three Josephson qubits setup, GHZ 
states of the form ^(|000) + e**|222)), g| can be obtained by applying two Q gates, the 
first one to b and a and the second one to b and c (with b always as the control qubit). 
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FIG. 1. (Top) Schematic representation of the proposed setup. The left junctions and the 
SQUID allow for the coherent manipulations to prepare the state. The right junctions provide the 
measurement setup. (Bottom) Procedure to prepare and detect the entanglement. The procedure 
in repeated as in the experiment in Ref 
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FIG. 2. Time Average of the current correlation function over the time interval At = 3r 1 as 
a function of the rotating angle 9 for the case of the mixture (left) and the singlet (right). The 



parameters in both figures are Ej 
E c h- 
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